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A digital equalizer for data transmission linearly combines a sequence 
of samples [b(r + kt), k = 0, ±1, • • -~\ofthe received data wave to miti- 
gate the effects of inter symbol interference and noise. A natural question is, 
How will the performance of such a system depend on the timing phase t, 
g t g T? 

We examine this problem in considerable detail for an infinite equalizer 
using a mean-square measure of performance. Excess bandwidth results 
in a significant difference in performance between the best and worst timing 
phase. Under practical noise conditions, it is estimated that the excess 
bandwidth must be down to 1 or 2 percent before the timing effect becomes 
insignificant. With a 10-percent roll-off, a 3-dB penalty can be incurred 
by choosing a bad timing epoch. 

Our main result is that under conditions likely to be encountered on 
channels similar to voiceband telephone channels, the optimum sampling 
instants will be accurately approximated by the consecutive maxima and 
minima of the sine wave that result when an alternating sequence of posi- 
tive and negative pulses {dotting sequence) is transmitted. There are no 
additional local minima of the minimum mean-square error as the timing 
phase is varied. 

I. INTRODUCTION 

We consider a noisy, random, baseband pulse train, 

b(t) = £ a„x{t - nT) + n(t), (1) 

where x (t) is the channel (or channel front-end filter) impulse response 
and the a„ are independent binary data which take values ±1 with 
equal probability. The additive zero mean gaussian noise process is 
denoted by n(t). Under very ideal conditions, the pulse x(t) would be 
of the form (sin rt/T)/(irt/T), and sampling (1) at time t = kT would 
yield the quantity a* + n(kT). Under more realistic conditions, the 
brick-wall shape of the above pulse is difficult to approximate in prac- 
tice, and a smoother characteristic in the frequency domain is taken as 
the ideal. Using some excess bandwidth (i.e., the Fourier transform 
extends beyond ir/T rad/s) still, in principle, allows the ideal set of 
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samples expressed by 

x(0) = 1, x(kT) = 0, k*0; 

but this ideal set is never exactly attained in practice because of un- 
known channel distortion. One possible measure of distortion could be 

E x 2 (kT) 

k*0 

X 2 (0) 

Actually, since there is now nothing special about choosing t = as 
the sampling instant for x(t), an even more appropriate measure 
would be 

E * 2 (r + kT) 
min*-" 



x\t) 

The above measure allows us to choose a best sampling epoch r before 
declaring how badly the signal has been distorted. However, any direct 
consideration of the above type of criterion seems to result in con- 
siderable mathematical difficulty. Also, in practice, one is more in- 
terested in the situation where the set of digital samples \x(t + kT)), 
or rather their noisy versions [b(r + kT)}, are linearly combined by 
an equalizer attempting to undo the channel distortion. It is the dis- 
tortion at the output of the equalizer that is of practical interest. 
Remarkably enough, it is this distortion measured at the output of an 
infinite equalizer that is most easily studied with regard to its properties 
as a function of timing phase. A discussion of this dependence follows. 
We begin by making more precise some of our general remarks. 

The received signal b(t) is to be equalized by passing it through a 
(2N + 1) tap transversal filter to yield for the signal portion of the 
output a waveform 

q(t) = E a n h(t - nT), (2) 

where 

h(t) - E c n x(t - nT). (3) 

-N 

The coefficients c„ are the tap weights of the equalizer and are adjusted 
to optimize some measure of performance. Since decisions are based 
on the sampled output q(t), one usually discusses the problem by fixing 
a sampling epoch r and denotes the sampled values of x(t) and h(t) 
by Xk and hk, respectively. That is, 

x(t + kT) m x k 

N k = ■• • -1,0, 1, •••. (4) 

h{ T + kT) = h k = E c n Xk-n 

n — N 
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We assume throughout that the tap weights are chosen to minimize 
the expression 

& 2 = f ' K + (h - l) 2 + <rl (5) 

where al is the sampled noise variance at the output of the equalizer, 
and the prime on the summation indicates deletion of the n = term. 
Letting R n = R-n be the (discrete) correlation function of the sampled 
noise at the input to the equalizer, we see that 

AT 

ol = L cjRj-kCk (6) 

i.k=-N 



and 



& = £ Ci(B/_* + A/-*)C* -2D CnX- n + 1, (7) 

j,k=-N n — N 



where the channel correlation matrix, 

CO 

A„- m = £ Xk-nXk-m, (8) 



k' 



has been introduced. The optimum tap settings (c opt ) n are determined 
by differentiating S 2 with respect to the tap gains. This yields a linear 
system of equations for the best tap gains which, in an obvious matrix- 
vector notation, reads 

{A + R)c opi = x. (9) 

Using this in (7), we then obtain for the minimum mean-square error 

SLn = 1 - ho. (10) 

Finally, we may state our problem. The above discussion of known 
results was in terms of a fixed timing epoch r. Here we shall determine 
for some special but important cases the value of t, which will minimize 
(Sfnm, as defined above. Instead of resorting to a direct numerical in- 
version of (9), we introduce the approximation of treating the equalizer 
as infinite. This approximation aids the analytical inversion of (9) and 
thereby directly leads to useful insights concerning the proper timing 
epoch. It will be seen that in the N = « limit, timing recovery is only a 
problem when the received pulse x(t) has a bandwidth greater than 
1/27* Hz, a result that could have been guessed from the sampling 
theorem. With the practical situation of excess bandwidth, there is the 
possibility that a badly chosen timing phase will considerably de- 
grade performance. Only when there is very little excess bandwidth, 
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say of the order of one percent, is the timing unimportant in practice. 
Our main result can be best stated in terms of the pure sine wave that 
appears at the channel output (sampler input) when the sequence 
a n = ( — 1)" (dotting sequence) is transmitted. The optimum sampling 
instants are, under certain conditions, the consecutive maxima and 
minima of this wave. Furthermore, if the timing phase is varied, no 
extraneous local minima of the minimum mean-square error (regarded 
as a function of timing phase) will appear. 

II. MATHEMATICAL DETERMINATION OF OPTIMUM TIMING PHASE 

On comparing (8) and (9), we see that the infinite matrix equation is 
in the nature of a convolution and, therefore, is best solved by trans- 
form methods. Thus, given a sequence \x„}, we define a transform 

X eq (w) = T f, Znexp (-iwnT), \w\ g |, (11) 

with the x k recoverable from X eq (w) in the obvious way : 

x k = ^[" X^(w)e*»"dw. (12) 

This formula shows that, from the impulse-response point of view, 
X e q(w) is, in fact, the Nyquist equivalent spectrum of x(t). That is, 
if X(w) is the Fourier transform of x(t), then 

X e<l (iv)= k i^x(w-k^, |w|g|. (13) 

In this transform language, the solution of (9) is* 

C( W ) ■ txim (14) 

In writing this, we have made use of (8) to write 

A(w) -ylX., (•)!■. (15) 

Equation (4) gives the relation 

H(w) = ±C(w)X eci (w) = |Xeq 1 ( ^ )|2 (16) 

fleqM + j, |-X"eq(w)| 2 



* We will not always use the subscript "eq" on the transform (11) when there does 
not exist an associated impulse response in a natural way. 
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and, hence, from (10) and the recovery formula for h , we have 



T r lT R ( w ) 



= L f 

27T J-(rlT 



) Bea(w) + ^ \X etl (w) 



dw. (17) 



Let us now explicitly put timing recovery into this solution for 8^. 
In evaluating the above formula, one would naturally calculate X(w), 
possibly by choosing t = as the peak of the pulse x(t), and calculate 
X eci (w) by (13). Any other timing epoch r can be obtained by replacing 
X(w) by X(w) exp (iwr). This corresponds to sampling the original 
x(t) at t = t instead of t = 0; in general, X eq (w) will depend on t. 

For further specifics, we limit our discussion to small but nonzero 
excess bandwidth, letting a denote the fraction of excess bandwidth. 
Thus, (1 + a)tr/T is the total bandwidth (in rad/s) occupied by the 
baseband pulse, a < 1. Further, we take the noise at the input to the 
equalizer to be independent from sample to sample and to have vari- 
ance a\\ thus, 

R eci (iv) = TV 2 , (18) 

Equation (17) thus becomes 

T r rlT dw 



2tT 7-(w/D _ 2 _|_ i 



T f'lTll-ai dw 



- a* i f 



D(l-o) 2 l Jl I Y" fmnW* 

+ «f f r^ as) 

^7T y + (,/r)(l-a) ! I i I V /,.A 12 

We note that the first term in the right-hand side of (19) is not affected 
by timing phase since the portion of \X e q(w)\ for |w>| ^ */T(l —a) 
is not affected, whereas the second term, on account of the foldover 
from excess bandwidth, is affected.* 

To obtain estimates regarding the magnitude of the effects of poor 
timing phase, we consider the case of an undisturbed pulse x(t) having 
a transform X(w) that is equal to T for ^ w ^ (1 - a)r/T and 
decreases linearly to zero for (1 — a)ir/T S w ^ (1 + a)r/T, with 
X(—w) = X(w). If this pulse is sampled at t = 0, the "main bang" 
has value unity, and there is no intersymbol interference. From (19), 
the minimum mean-square error is <r 2 t (l + a*) £zi a\ for small noise. 



* Of course, if there is no noise, the equalizer inverts X eQ (w) (if the latter is not zero) 
without any penalty, and the question of timing is moot. 
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If, however, it is sampled at t = T/2, we have 



T\ < w ^ (1 - a) | 

(20) 






and, hence, with the infinite equalizer doing the best it can, we obtain 
from (19) for small noise 

SLn = (1 - *)o\ + a<ri tan" 1 - (21) 

W(l -a)^+|a<r,. 

If for binary transmission we assume a 10 or 15 percent roll-off (a = 0.1 
or 0.15) and a\ = 0.01 to 0.04, we see that a degradation of from 2 dB 
to 4 dB might be encountered. Only if the roll-off is of the order of 
1 or 2 percent is the effect of timing negligible for the value of noise 
power considered here. 

We now concentrate further attention on the second term of (19). 
We introduce a variable 

>-=[•-*] (22) 

for ir/T(l - a) ^ w ^ t/T(1 + <*); thus, -1 £ p £ 1. When w and 
n are related by (22), we define a nonnegative amplitude function R(w) 
and a real phase function <p(n) by the relation 

X(w) = BGb)«OW, B(m) ^ 0. (23) 

Further, we introduce the even and odd parts of amplitude and phase 
about the Nyquist frequency w = t/T(ji = 0) : 



R(n) = BM + Ro(n) 
RM m B.(- M );ieo(M) = -Ao(-m); 

<p(v) = V«(m) + ^o(m) 

ipo(ti) = (p t {— n);<p (fj.) = —tpo(-n). 



(24) 
(25) 



Thus, for — 1 5g /j ^ 0, we calculate 

|X eq (/x) I 2 = 4[(# 2 - fl 2 ) cos 2 a + BQ (26) 

where, for notational simplicity, the ji dependence of the function on 
the right has not been shown. Using (26), the second term S of (19) is 
written 



S = 



atr t 



p 1- (27) 

-/-l 2 i ^ r/r>2 ri2\ •> i r>2~i 



* 2 < + f-2 U R * ~ & COs2 * + fl ^ 
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A number of remarks concerning (27) are in order. 

(i) The odd component of phase distortion about the Nyquist 
frequency does not influence "foldover" or the timing recovery 
problem. 
(ii) From the positivity of R(ji), we note that (R 2 e - R$) ^ 0. We 
assume that [# 2 (m) — fljJGO] > ° on some interval contained 
in - 1 ^ fi ^ 0. 
(Hi) Relative to an arbitrarily selected sampling instant, called zero, 
the timing phase r enters <p e through 

V» = y r + £.(/*), (28) 

where £„(m) is the phase characteristic for t = 0. It is evident 
that the mean-square error is a periodic function of r with 
period T. Physically, if r is displaced by an integer (k) mul- 
tiple of T, the "center tap" will be shifted k units down the 
equalizer. 
(iv) In many situations, ip e can, to a first approximation, be treated 
as a constant. A recent survey 1 shows that an "average" tele- 
phone channel has an envelope delay that changes by 312 /xs 
from 2850 Hz to 3000 Hz. Assuming this to be an excess band- 
width region and the change to be due entirely to a quadratic 
term in the phase curve yields a phase change of 8.4 degrees.* 

We now come to the main point of this section, which is to recognize 
where the optimum timing phase is. This will be done with the sim- 
plifying assumption 

<p e = 8 = const. (29) 

Since S is a smooth periodic function of t, we look for a minimum by 
examining solutions of 

^ = 0, (30) 

dr 

that is, 

sin 2<p e / -f j n - U- (61) 

J ~ 1 Wi + fi [ («! - B8) cos 2 <p e + R%] 

From (ii) we deduce that the integral is strictly positive and, hence, 
the only solutions are those of 

sin2^ e = 0. (32) 



* This is meant to be a rough estimate only. No attempt was made to separate out 
even and odd contributions to the change in phase. Low-end envelope delay is more 
severe and this separation should be carried out so as to obtain a not-too-pessimistic 
answer. 
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Thus, because of periodicity the only solutions t of interest may be 
written 

| r + 6 = (33a) 



or 



| r + + \ = 0. (33b) 

By taking second derivatives of S with respect to t, we see that the 
first of these, namely (33a), corresponds to the minimum of S while 
the second solution is a maximum. In addition to giving a description 
of optimum timing phase, our discussion has also shown that for 
<p e = const, there are no other local minima of <§mm as t is varied. * Thus, 
a gradient search for the optimum r is feasible. More precisely, we 
have shown this for the quadratic Nyquist I criterion (5), and have 
assumed that the search is slow in the sense that the equalizer "settles 
down" before the next step in t is taken. 

Equation (33a) provides a mathematical description of the optimum 
timing phase t in terms of the phase of the channel at the Nyquist 
frequency t/T. It is useful to give a more physical interpretation of 
this result. Consider the transmission of the special sequence a„ 
= ( — 1)". From (1), we receive the signal 

t (-l)»x(t - nT), 

which has as its Fourier transform 

X(w) £ (-l)"exp (-iwnT) 

= X(w) f; exp (-iw2kT) - X{w) exp (-iwT) £ exp (-iw2kT) 
= Y x ( w )tj( w -Y)-Y X{w)exp( - iwT) tj( w -Y) 

Hence, if X(w) is band-limited to a bandwidth less than 3tr/T, and if 
we write X(t/L) = A exp (id), the received signal is 



4 , 
j, A cos 



Comparing this with (32) and (33a), we see that the optimum timing 



* See the appendix for a more general condition under which the existence of a 
unique minimum is guaranteed. 
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phase is at the successive maxima and minima of the tone that results 
when a dotting sequence is transmitted. 

Finally, we note that corrections for a small, varying, even compo- 
nent of phase distortion about the Nyquist frequency are easy to make. 
If we define 7 by 

£t + = t. (34) 

where, again, is the phase at the Nyquist frequency, and if we write 
[recalling (28)], 

<pM = 7 + «(/*), e(0) = 0, 



where c(m) is assumed small, then 

(35) 






i: 



(m-RD dM 



a* + (4/T*)R 2 e \> 



where 7* is the approximate value of 7 for optimum timing phase t. 
We see 7* is bounded in magnitude by max e(/x), but in general it will 
be a fraction of this. Assuming max e(/x) corresponds to 9 degrees, this 
rough bound suggests that a correction of at most 5 percent of a sym- 
bol interval would be required. 

For a more refined estimate we take the following 

R e (n) = 1, R (n) = M, and e(/x) = cm 2 - (36) 

Note that R\ - Rl = {R. - Ro)(R. + Ro) = R(ji)R(-v). Since R(l) 
will vanish by the design of the signal, R 2 e — Rl will be zero at n = — 1 
and, thus, no multiplicative constant is included in ftoU)- Using (36) 
in (35) yields 



J (1 - mV<*m 

£ (1 - M 2 )dM 



7*«-e^- - " -i' (37) 



Thus an estimate of \ of the even part of the phase variation over the 
roll-off band results for the correction term to (34). This would appear 
to be a negligible effect. 

III. DISCUSSION FOR PASSBAND TRANSMISSION 

We now generalize the foregoing discussion to treat a qam passband 
signal format. We assume a transmitted signal 

[L a n p{t - nT)2 cos w c t - [£ b n p(t - nT)~] sin w c t, (38) 
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where (a„) and (&„) are independent digits. Let the impulse response 
of the linear distortion medium m (t) be written as 

m(t) = 2Fi(0 cos w e t - 2F 2 (t) sin w e t. (39) 

Further, denote the convolution of p(t) with F\(t) and F»(t) by x(t) 
and y(0, i.e., 

p(t)*Fi(t) = x(t) p(t)*F 2 (t) = y{t). (40) 

Then the received signal is 

r(t) = cos w c t[_Y, a n x{t - nT) - £ 6„2/(< - nT 1 )] 

- sin »■*[£ a„y(* - nT) + £ 6„x(t - nT)] + n(0, (41) 

where n{t) represents additive gaussian noise. The waveform r(t) is then 
linearly filtered by a passband equalizer, whose output can be de- 
scribed by replacing x(t) by g(t) and y(t) by h(t), where 

0(0 = £ c n x(t - nT) - £ d»y(l - nT) 

-AT -AT 

and (42) 

MO = £ d n x(t -nT) + £ c„jK« - nr). 

-AT -AT 

Further, if the input noise n(t) has correlation R(t), then the output 
noise has variance 

<rg = c+flc + d+ita, (43) 

where Ra = R[(i — j)T~\. We adopt as our criterion 

& - £ 2 , + (1 - <7o) 2 + £ /i? + o-B (44) 

and seek to minimize this quantity. Analogously with the baseband 
case we define two baseband autocorrelation functions {^corresponding 
to x(t) and y(t) in (41)] by 

00 

**-nm «■« — m 2-t * j—nXj—m 

j =—oo 

and (45) 

00 

** nm = Jjn—m = 2^ V J—nU j—m< 

j— 00 

Further, a cross-correlation matrix between x(t) and y(t), namely, 

Knm = Kn-m = —K mn = £ (x j-„y j- m — Xj-. m yj- n ). (46) 

i 
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If, as before, we write (x) n = x- n , (y)n = y~ n , then 

& = c+(A + B)c + d+(A + B)d - 2c + Kd - 2c-x + 2d-y 

+ 1 + c+Rc + d+Rd. (47) 

If one defines the matrix Q = A -\- B + R, the equations for the 
optimum taps are 

£ (Qn- m C m - K n - m d m ) = X_ n 



E (Qn_ m d m + K n - m c m ) = -y-„ 
and 



(48) 



£Lm = 1 - (go) opt = 1 - c 0P t -a; + dopt-y. (49) 

We introduce the same transform as before, noting 

Q(w) = |Xeq ^ )|2 + |Feq ^ )|2 + R(w) = Q{-w) (50) 



and 



K(w) = \ lXl,(w)Y eci (w) - Ze^^n,^)] = -K(-w). (51) 



Straightforward solution of (48) and (49) yields 

T ni* Q(w)R(w) 
6m,n _ 2* j_ (wr) Q*(w) + tf 2 (i/;) aw 

_T /" /7, [QM - ig(u;)]fl(«;) , 
" 2t J_ {T/ r) Q 2 (™) + K*(w) 

2tt J- ( ,/d Q(w) + ia(w) 

where the oddness of K(w) [and evenness of Q(w) and R{w)~\ has been 
used. If we formally introduce S(w) by 

S(w) = X(w) + iY(w), (53) 

and 

S eQ » = *ea» + <7«(tf)i (54) 

we see that 

!«•*(») I" = LXa(«0 + *T^(»)XJEi(to) - iHt(tp)] 

= |X eq (w;)| 2 + |r ea »| 2 

+ »F^(w)JR,(w) - iX^(w)Yl q (w) 

= T[A(u>) +BW + ii? (to)} (55) 

or, briefly, 

Q + iK = i^ + R. (56) 
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Thus, 

02. T f" T R ( w ) J /e*\ 

5Ln = o- / ■» .v., dw. (57) 

This formula has a close analogy to the baseband formula when we 
realize that the overall channel (transmitting and distorting filter) 
has frequency characteristic 

C(w) = hS(w - w c ) + $S*(-w - w c ). (58) 

Thus, in terms of the overall channel, (57) has a straightforward 
interpretation. One uses a formula similar to the baseband version (17) 
except that now one forms an equivalent Nyquist characteristic using 
twice the positive frequency part of the passband channel and treats 
the carrier frequency as zero. It is essential to note that the transfer 
characteristic from which S eci (w) is formed no longer needs to have 
even amplitude and odd phase about "zero" as would be required in 
the baseband case. In particular, the envelope delays at the two band 
edges can have an effect on timing phase. In practical cases, the carrier 
frequency can be chosen to give equal envelope delays, thus minimizing 
the effects of odd phase about the band edges. However, in doing this, 
the amplitude effects may have worsened the situation. The main 
point to be emphasized here is the interplay of carrier placement and 
bit timing with an excess bandwidth system. 

Regarding the description of optimum timing phase, let #* be the 
channel phase at w = w c ± tr/T, and assume these phases to be con- 
stant over their respective roll-off regions. Then, for any reasonable 
amplitude distortions, the unique best sampling instant t is given by 

§r + *(*+- <p~) =0. (59) 

Under these conditions, we again have that there exists but one local 
minimum. 
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APPENDIX 

Uniqueness of Local Minimum of Mean-Square Error for Small ip.(x) 

The interesting portion of the mean-square error is given by (27), 
which is of the form 
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g{y) = /, a(x) + b(x) cos' [7 + *»(x)J ' (60) 

with o(«) and b(x) positive. If <p(x) is always small, then any local 
minimum of (60) has to occur when 7 is small. Thus, if | <p(x) \ < x/8, 
we see that the derivative 

K ^ T b(x) sin [2 7 + 2^(a;)] (61) 

g (7) = Jo {a(x)+b(x)co S *ly+<p(x)]\> ^ 

can only vanish if 1 7 1 < r/B or | tt - 7 1 < */8, for otherwise the 
integrand always has the same sign. Let us investigate the solution 
in 1 7 1 < r/8 in more detail. Taking another derivative gives 

fit \ - 9 f 1 A, b C0S (27 + 2<p) 
g (7) -2j Q <fa [o + 6 COB i ( 7 + rfj 

; [a + 6 cos 2 (7 + p)] 3 

which is positive. Hence, there can be only one solution of g'{y) = 
for I7I < ir/S (we cannot have two minima without a maximum in 
between). Likewise, if | <p(x) | < t/8, there is precisely one maximum 
of g(y), and it occurs in the range \v — y\ < r/8. 
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